Deciphering the spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer.

Although heterogeneity of FAP+ Cancer-Associated Fibroblasts (CAF) has been described in breast cancer, their plasticity and spatial distribution remain poorly understood. Here, we analyze trajectory inference, deconvolute spatial transcriptomics at single-cell level and perform functional assays to generate a high-resolution integrated map of breast cancer (BC), with a focus on inflammatory and myofibroblastic (iCAF/myCAF) FAP+ CAF clusters. We identify 10 spatially-organized FAP+ CAF-related cellular niches, called EcoCellTypes, which are differentially localized within tumors. Consistent with their spatial organization, cancer cells drive the transition of detoxification-associated iCAF (Detox-iCAF) towards immunosuppressive extracellular matrix (ECM)-producing myCAF (ECM-myCAF) via a DPP4- and YAP-dependent mechanism. In turn, ECM-myCAF polarize TREM2+ macrophages, regulatory NK and T cells to induce immunosuppressive EcoCellTypes, while Detox-iCAF are associated with FOLR2+ macrophages in an immuno-protective EcoCellType. FAP+ CAF subpopulations accumulate differently according to the invasive BC status and predict invasive recurrence of ductal carcinoma in situ (DCIS), which could help in identifying low-risk DCIS patients eligible for therapeutic de-escalation.

Residual ANTXR1+ myofibroblasts after chemotherapy inhibit anti-tumor immunity via YAP1 signaling pathway.

Although cancer-associated fibroblast (CAF) heterogeneity is well-established, the impact of chemotherapy on CAF populations remains poorly understood. Here we address this question in high-grade serous ovarian cancer (HGSOC), in which we previously identified 4 CAF populations. While the global content in stroma increases in HGSOC after chemotherapy, the proportion of FAP+ CAF (also called CAF-S1) decreases. Still, maintenance of high residual CAF-S1 content after chemotherapy is associated with reduced CD8+ T lymphocyte density and poor patient prognosis, emphasizing the importance of CAF-S1 reduction upon treatment. Single cell analysis, spatial transcriptomics and immunohistochemistry reveal that the content in the ECM-producing ANTXR1+ CAF-S1 cluster (ECM-myCAF) is the most affected by chemotherapy. Moreover, functional assays demonstrate that ECM-myCAF isolated from HGSOC reduce CD8+ T-cell cytotoxicity through a Yes Associated Protein 1 (YAP1)-dependent mechanism. Thus, efficient inhibition after treatment of YAP1-signaling pathway in the ECM-myCAF cluster could enhance CD8+ T-cell cytotoxicity. Altogether, these data pave the way for therapy targeting YAP1 in ECM-myCAF in HGSOC.

Molecular underpinnings and environmental drivers of loss of heterozygosity in Drosophila intestinal stem cells.

During development and aging, genome mutation leading to loss of heterozygosity (LOH) can uncover recessive phenotypes within tissue compartments. This phenomenon occurs in normal human tissues and is prevalent in pathological genetic conditions and cancers. While studies in yeast have defined DNA repair mechanisms that can promote LOH, the predominant pathways and environmental triggers in somatic tissues of multicellular organisms are not well understood. Here, we investigate mechanisms underlying LOH in intestinal stem cells in Drosophila. Infection with the pathogenic bacteria, Erwinia carotovora carotovora 15, but not Pseudomonas entomophila, increases LOH frequency. Using whole genome sequencing of somatic LOH events, we demonstrate that they arise primarily via mitotic recombination. Molecular features and genetic evidence argue against a break-induced replication mechanism and instead support cross-over via double Holliday junction-based repair. This study provides a mechanistic understanding of mitotic recombination, an important mediator of LOH, and its effects on stem cells in vivo.

KDM6B drives epigenetic reprogramming associated with lymphoid stromal cell early commitment and immune properties.

Mature lymphoid stromal cells (LSCs) are key organizers of immune responses within secondary lymphoid organs. Similarly, inflammation-driven tertiary lymphoid structures depend on immunofibroblasts producing lymphoid cytokines and chemokines. Recent studies have explored the origin and heterogeneity of LSC/immunofibroblasts, yet the molecular and epigenetic mechanisms involved in their commitment are still unknown. This study explored the transcriptomic and epigenetic reprogramming underlying LSC/immunofibroblast commitment. We identified the induction of lysine demethylase 6B (KDM6B) as the primary epigenetic driver of early immunofibroblast differentiation. In addition, we observed an enrichment for KDM6B gene signature in murine inflammatory fibroblasts and pathogenic stroma of patients with autoimmune diseases. Last, KDM6B was required for the acquisition of LSC/immunofibroblast functional properties, including the up-regulation of CCL2 and the resulting recruitment of monocytes. Overall, our results reveal epigenetic mechanisms that participate in the early commitment and immune properties of immunofibroblasts and support the use of epigenetic modifiers as fibroblast-targeting strategies in chronic inflammation.

Imaging and multi-omics datasets converge to define different neural progenitor origins for ATRT-SHH subgroups.

Atypical teratoid rhabdoid tumors (ATRT) are divided into MYC, TYR and SHH subgroups, suggesting diverse lineages of origin. Here, we investigate the imaging of human ATRT at diagnosis and the precise anatomic origin of brain tumors in the Rosa26-CreERT2::Smarcb1flox/flox model. This cross-species analysis points to an extra-cerebral origin for MYC tumors. Additionally, we clearly distinguish SHH ATRT emerging from the cerebellar anterior lobe (CAL) from those emerging from the basal ganglia (BG) and intra-ventricular (IV) regions. Molecular characteristics point to the midbrain-hindbrain boundary as the origin of CAL SHH ATRT, and to the ganglionic eminence as the origin of BG/IV SHH ATRT. Single-cell RNA sequencing on SHH ATRT supports these hypotheses. Trajectory analyses suggest that SMARCB1 loss induces a de-differentiation process mediated by repressors of the neuronal program such as REST, ID and the NOTCH pathway.

Multi-omics comparison of malignant and normal uveal melanocytes reveals molecular features of uveal melanoma.

Uveal melanoma (UM) is a rare cancer resulting from the transformation of melanocytes in the uveal tract. Integrative analysis has identified four molecular and clinical subsets of UM. To improve our molecular understanding of UM, we performed extensive multi-omics characterization comparing two aggressive UM patient-derived xenograft models with normal choroidal melanocytes, including DNA optical mapping, specific histone modifications, and DNA topology analysis using Hi-C. Our gene expression and cytogenetic analyses suggest that genomic instability is a hallmark of UM. We also identified a recurrent deletion in the BAP1 promoter resulting in loss of expression and associated with high risk of metastases in UM patients. Hi-C revealed chromatin topology changes associated with the upregulation of PRAME, an independent prognostic biomarker in UM, and a potential therapeutic target. Our findings illustrate how multi-omics approaches can improve our understanding of tumorigenesis and reveal two distinct mechanisms of gene expression dysregulation in UM.

The mutational landscape of skull base and spinal chordomas and the identification of potential prognostic and theranostic biomarkers.

OBJECTIVE: Chordomas are rare bone neoplasms characterized by a high recurrence rate and no benefit from any approved medical treatment to date. However, the investigation of molecular alterations in chordomas could be essential to prognosticate, guide clinical decision-making, and identify theranostic biomarkers. The aim of this study was to provide a detailed genomic landscape of a homogeneous series of 64 chordoma samples, revealing driver events, theranostic markers, and outcome-related genomic features. METHODS: The authors conducted whole-exome sequencing (WES), targeted next-generation sequencing, and RNA sequencing of 64 skull base and spinal chordoma samples collected between December 2006 and September 2020. Clinical, histological, and radiological data were retrospectively analyzed and correlated to genetic findings. RESULTS: The authors identified homozygous deletions of CDKN2A/2B, PIK3CA mutations, and alterations affecting genes of SWI/SNF chromatin remodeling complexes (PBRM1 and ARID1A) as potential theranostic biomarkers. Using matched germline WES, they observed a higher frequency of a common genetic variant (rs2305089; p.(Gly177Asp)) in TBXT (97.8%, p < 0.001) compared to its distribution in the general population. PIK3CA mutation was identified as an independent biomarker of short progression-free survival (HR 10.68, p = 0.0008). Loss of CDKN2A/2B was more frequently observed in spinal tumors and recurrent tumors. CONCLUSIONS: In the current study, the authors identified driver events such as PBRM1 and PIK3CA mutations, TBXT alterations, or homozygous deletions of CDKN2A/2B, which could, for some, be considered potential theranostic markers and could allow for identifying novel therapeutic approaches. With the aim of a future biomolecular prognostication classification, alterations affecting PIK3CA and CDKN2A/2B could be considered as poor prognostic biomarkers.

Epigenetically controlled tumor antigens derived from splice junctions between exons and transposable elements.

Oncogenesis often implicates epigenetic alterations, including derepression of transposable elements (TEs) and defects in alternative splicing. Here, we explore the possibility that noncanonical splice junctions between exons and TEs represent a source of tumor-specific antigens. We show that mouse normal tissues and tumor cell lines express wide but distinct ranges of mRNA junctions between exons and TEs, some of which are tumor specific. Immunopeptidome analyses in tumor cell lines identified peptides derived from exon-TE splicing junctions associated to MHC-I molecules. Exon-TE junction-derived peptides were immunogenic in tumor-bearing mice. Both prophylactic and therapeutic vaccinations with junction-derived peptides delayed tumor growth in vivo. Inactivation of the TE-silencing histone 3-lysine 9 methyltransferase Setdb1 caused overexpression of new immunogenic junctions in tumor cells. Our results identify exon-TE splicing junctions as epigenetically controlled, immunogenic, and protective tumor antigens in mice, opening possibilities for tumor targeting and vaccination in patients with cancer.

Single-cell transcriptomics reveals shared immunosuppressive landscapes of mouse and human neuroblastoma.

BACKGROUND: High-risk neuroblastoma is a pediatric cancer with still a dismal prognosis, despite multimodal and intensive therapies. Tumor microenvironment represents a key component of the tumor ecosystem the complexity of which has to be accurately understood to define selective targeting opportunities, including immune-based therapies. METHODS: We combined various approaches including single-cell transcriptomics to dissect the tumor microenvironment of both a transgenic mouse neuroblastoma model and a cohort of 10 biopsies from neuroblastoma patients, either at diagnosis or at relapse. Features of related cells were validated by multicolor flow cytometry and functional assays. RESULTS: We show that the immune microenvironment of MYCN-driven mouse neuroblastoma is characterized by a low content of T cells, several phenotypes of macrophages and a population of cells expressing signatures of myeloid-derived suppressor cells (MDSCs) that are molecularly distinct from the various macrophage subsets. We document two cancer-associated fibroblasts (CAFs) subsets, one of which corresponding to CAF-S1, known to have immunosuppressive functions. Our data unravel a complex content in myeloid cells in patient tumors and further document a striking correspondence of the microenvironment populations between both mouse and human tumors. We show that mouse intratumor T cells exhibit increased expression of inhibitory receptors at the protein level. Consistently, T cells from patients are characterized by features of exhaustion, expressing inhibitory receptors and showing low expression of effector cytokines. We further functionally demonstrate that MDSCs isolated from mouse neuroblastoma have immunosuppressive properties, impairing the proliferation of T lymphocytes. CONCLUSIONS: Our study demonstrates that neuroblastoma tumors have an immunocompromised microenvironment characterized by dysfunctional T cells and accumulation of immunosuppressive cells. Our work provides a new and precious data resource to better understand the neuroblastoma ecosystem and suggest novel therapeutic strategies, targeting both tumor cells and components of the microenvironment.

Enrichment of centromeric DNA from human cells.

Centromeres are key elements for chromosome segregation. Canonical centromeres are built over long-stretches of tandem repetitive arrays. Despite being quite abundant compared to other loci, centromere sequences overall still represent only 2 to 5% of the human genome, therefore studying their genetic and epigenetic features is a major challenge. Furthermore, sequencing of centromeric regions requires high coverage to fully analyze length and sequence variations, and this can be extremely costly. To bypass these issues, we have developed a technique, named CenRICH, to enrich for centromeric DNA from human cells based on selective restriction digestion and size fractionation. Combining restriction enzymes cutting at high frequency throughout the genome, except within most human centromeres, with size-selection of fragments >20 kb, resulted in over 25-fold enrichment in centromeric DNA. High-throughput sequencing revealed that up to 60% of the DNA in the enriched samples is made of centromeric repeats. We show that this method can be used in combination with long-read sequencing to investigate the DNA methylation status of certain centromeres and, with a specific enzyme combination, also of their surrounding regions (mainly HSATII). Finally, we show that CenRICH facilitates single-molecule analysis of replicating centromeric fibers by DNA combing. This approach has great potential for making sequencing of centromeric DNA more affordable and efficient and for single DNA molecule studies.

Dichotomy in Neutralizing Antibody Induction to Peptide-Conjugated Vaccine in Squalene Emulsion Contrast With Aluminum Hydroxide Formulation.

W614A-3S peptide is a modified 3S motif of the HIV-gp41 (mutation W614A). We previously detected the presence of natural neutralizing antibodies directed against W614A-3S peptide (NAbs) in long-term non-progressor HIV+ patients. Here, we compared the efficacy of W614A-3S peptide formulated in either squalene emulsion (SQE) or in aluminum hydroxide (Alum) in inducing broadly-NAbs (bNAbs). Rabbit and mouse models were used to screen the induction of bNAbs following 4 immunizations. SQE was more efficient than Alum formulation in inducing W614A-3S-specific bNAbs with up to 67%-93% of HIV strains neutralized. We then analyzed the quality of peptide-specific murine B cells by single-cell gene expression by quantitative reverse transcription-PCR and single-cell V(D)J sequencing. We found more frequent germinal center B cells in SQE than in Alum, albeit with a different gene expression profile. The V(D)J sequencing of W614A-3S-specific BCR showed significant differences in BCR sequences and validates the dichotomy between adjuvant formulations. All sixteen BCR sequences which were cloned were specific of peptide. Adjuvant formulations of W614A-3S-peptide-conjugated immunogen impact the quantity and quality of B cell immune responses at both the gene expression level and BCR sequence.

PD-L1 and ICOSL discriminate human Secretory and Helper dendritic cells in cancer, allergy and autoimmunity.

Dendritic cells (DC) are traditionally classified according to their ontogeny and their ability to induce T cell response to antigens, however, the phenotypic and functional state of these cells in cancer does not necessarily align to the conventional categories. Here we show, by using 16 different stimuli in vitro that activated DCs in human blood are phenotypically and functionally dichotomous, and pure cultures of type 2 conventional dendritic cells acquire these states (termed Secretory and Helper) upon appropriate stimuli. PD-L1highICOSLlow Secretory DCs produce large amounts of inflammatory cytokines and chemokines but induce very low levels of T helper (Th) cytokines following co-culturing with T cells. Conversely, PD-L1lowICOSLhigh Helper DCs produce low levels of secreted factors but induce high levels and a broad range of Th cytokines. Secretory DCs bear a single-cell transcriptomic signature indicative of mature migratory LAMP3+ DCs associated with cancer and inflammation. Secretory DCs are linked to good prognosis in head and neck squamous cell carcinoma, and to response to checkpoint blockade in Melanoma. Hence, the functional dichotomy of DCs we describe has both fundamental and translational implications in inflammation and immunotherapy.

H3K27me3 conditions chemotolerance in triple-negative breast cancer.

The persistence of cancer cells resistant to therapy remains a major clinical challenge. In triple-negative breast cancer, resistance to chemotherapy results in the highest recurrence risk among breast cancer subtypes. The drug-tolerant state seems largely defined by nongenetic features, but the underlying mechanisms are poorly understood. Here, by monitoring epigenomes, transcriptomes and lineages with single-cell resolution, we show that the repressive histone mark H3K27me3 (trimethylation of histone H3 at lysine 27) regulates cell fate at the onset of chemotherapy. We report that a persister expression program is primed with both H3K4me3 (trimethylation of histone H3 at lysine 4) and H3K27me3 in unchallenged cells, with H3K27me3 being the lock to its transcriptional activation. We further demonstrate that depleting H3K27me3 enhances the potential of cancer cells to tolerate chemotherapy. Conversely, preventing H3K27me3 demethylation simultaneously to chemotherapy inhibits the transition to a drug-tolerant state, and delays tumor recurrence in vivo. Our results highlight how chromatin landscapes shape the potential of cancer cells to respond to initial therapy.

Kronos scRT: a uniform framework for single-cell replication timing analysis.

Mammalian genomes are replicated in a cell type-specific order and in coordination with transcription and chromatin organization. Currently, single-cell replication studies require individual processing of sorted cells, yielding a limited number (<100) of cells. Here, we develop Kronos scRT, a software for single-cell Replication Timing (scRT) analysis. Kronos scRT does not require a specific platform or cell sorting, which allows investigating large datasets obtained from asynchronous cells. By applying our tool to published data as well as droplet-based single-cell whole-genome sequencing data generated in this study, we exploit scRT from thousands of cells for different mouse and human cell lines. Our results demonstrate that although genomic regions are frequently replicated around their population average RT, replication can occur stochastically throughout S phase. Altogether, Kronos scRT allows fast and comprehensive investigations of the RT programme at the single-cell resolution for both homogeneous and heterogeneous cell populations.

Tissue-resident FOLR2+ macrophages associate with CD8+ T cell infiltration in human breast cancer.

Macrophage infiltration is a hallmark of solid cancers, and overall macrophage infiltration correlates with lower patient survival and resistance to therapy. Tumor-associated macrophages, however, are phenotypically and functionally heterogeneous. Specific subsets of tumor-associated macrophage might be endowed with distinct roles on cancer progression and antitumor immunity. Here, we identify a discrete population of FOLR2+ tissue-resident macrophages in healthy mammary gland and breast cancer primary tumors. FOLR2+ macrophages localize in perivascular areas in the tumor stroma, where they interact with CD8+ T cells. FOLR2+ macrophages efficiently prime effector CD8+ T cells ex vivo. The density of FOLR2+ macrophages in tumors positively correlates with better patient survival. This study highlights specific roles for tumor-associated macrophage subsets and paves the way for subset-targeted therapeutic interventions in macrophages-based cancer therapies.

PLK1 inhibition exhibits strong anti-tumoral activity in CCND1-driven breast cancer metastases with acquired palbociclib resistance.

A significant proportion of patients with oestrogen receptor (ER) positive breast cancers (BC) develop resistance to endocrine treatments (ET) and relapse with metastatic disease. Here we perform whole exome sequencing and gene expression analysis of matched primary breast tumours and bone metastasis-derived patient-derived xenografts (PDX). Transcriptomic analyses reveal enrichment of the G2/M checkpoint and up-regulation of Polo-like kinase 1 (PLK1) in PDX. PLK1 inhibition results in tumour shrinkage in highly proliferating CCND1-driven PDX, including different RB-positive PDX with acquired palbociclib resistance. Mechanistic studies in endocrine resistant cell lines, suggest an ER-independent function of PLK1 in regulating cell proliferation. Finally, in two independent clinical cohorts of ER positive BC, we find a strong association between high expression of PLK1 and a shorter metastases-free survival and poor response to anastrozole. In conclusion, our findings support clinical development of PLK1 inhibitors in patients with advanced CCND1-driven BC, including patients progressing on palbociclib treatment.

Transcriptional Programs Define Intratumoral Heterogeneity of Ewing Sarcoma at Single-Cell Resolution.

EWSR1-FLI1, the chimeric oncogene specific for Ewing sarcoma (EwS), induces a cascade of signaling events leading to cell transformation. However, it remains elusive how genetically homogeneous EwS cells can drive the heterogeneity of transcriptional programs. Here, we combine independent component analysis of single-cell RNA sequencing data from diverse cell types and model systems with time-resolved mapping of EWSR1-FLI1 binding sites and of open chromatin regions to characterize dynamic cellular processes associated with EWSR1-FLI1 activity. We thus define an exquisitely specific and direct enhancer-driven EWSR1-FLI1 program. In EwS tumors, cell proliferation and strong oxidative phosphorylation metabolism are associated with a well-defined range of EWSR1-FLI1 activity. In contrast, a subpopulation of cells from below and above the intermediary EWSR1-FLI1 activity is characterized by increased hypoxia. Overall, our study reveals sources of intratumoral heterogeneity within EwS tumors.

FLASH Irradiation Spares Lung Progenitor Cells and Limits the Incidence of Radio-induced Senescence.

PURPOSE: One of the main limitations to anticancer radiotherapy lies in irreversible damage to healthy tissues located within the radiation field. "FLASH" irradiation at very high dose-rate is a new treatment modality that has been reported to specifically spare normal tissue from late radiation-induced toxicity in animal models and therefore could be a promising strategy to reduce treatment toxicity. EXPERIMENTAL DESIGN: Lung responses to FLASH irradiation were investigated by qPCR, single-cell RNA sequencing (sc-RNA-Seq), and histologic methods during the acute wound healing phase as well as at late stages using C57BL/6J wild-type and Terc-/- mice exposed to bilateral thorax irradiation as well as human lung cells grown in vitro. RESULTS: In vitro studies gave evidence of a reduced level of DNA damage and induced lethality at the advantage of FLASH. In mouse lung, sc-RNA-seq and the monitoring of proliferating cells revealed that FLASH minimized the induction of proinflammatory genes and reduced the proliferation of progenitor cells after injury. At late stages, FLASH-irradiated lungs presented less persistent DNA damage and senescent cells than after CONV exposure, suggesting a higher potential for lung regeneration with FLASH. Consistent with this hypothesis, the beneficial effect of FLASH was lost in Terc-/- mice harboring critically short telomeres and lack of telomerase activity. CONCLUSIONS: The results suggest that, compared with conventional radiotherapy, FLASH minimizes DNA damage in normal cells, spares lung progenitor cells from excessive damage, and reduces the risk of replicative senescence.

Discovery of a new predominant cytosine DNA modification that is linked to gene expression in malaria parasites.

DNA cytosine modifications are key epigenetic regulators of cellular processes in mammalian cells, with their misregulation leading to varied disease states. In the human malaria parasite Plasmodium falciparum, a unicellular eukaryotic pathogen, little is known about the predominant cytosine modifications, cytosine methylation (5mC) and hydroxymethylation (5hmC). Here, we report the first identification of a hydroxymethylcytosine-like (5hmC-like) modification in P. falciparum asexual blood stages using a suite of biochemical methods. In contrast to mammalian cells, we report 5hmC-like levels in the P. falciparum genome of 0.2-0.4%, which are significantly higher than the methylated cytosine (mC) levels of 0.01-0.05%. Immunoprecipitation of hydroxymethylated DNA followed by next generation sequencing (hmeDIP-seq) revealed that 5hmC-like modifications are enriched in gene bodies with minimal dynamic changes during asexual development. Moreover, levels of the 5hmC-like base in gene bodies positively correlated to transcript levels, with more than 2000 genes stably marked with this modification throughout asexual development. Our work highlights the existence of a new predominant cytosine DNA modification pathway in P. falciparum and opens up exciting avenues for gene regulation research and the development of antimalarials.